o 
o 
o 

(N 

a* 

(N 



A BROAD Fe Ka EMISSION LINE IN THE X-RAY 
SPECTRUM OF THE QUASAR 3C 273 

Tahir Yaqoob 1,2 , Peter Serlemitsos 1 



Accepted for Publication in the Astrophysical Journal Letters (9/25/00; submitted 

8/9/2000) 



> 

in 
m 
^r- 

o\ 

O ! Received ; accepted 

O 
O 

a: 
6 

in 
c3 



% 



laboratory for High Energy Astrophysics, NASA/Goddard Space Flight Center, 
Greenbelt, MD 20771, USA. 

2 Joint Center for Astrophysics, Physics Department, University of Maryland, Baltimore 
County, Baltimore, 1000 Hilltop Circle, MD21250. 



-2- 



ABSTRACT 



The Fe Ka line, an important physical diagnostic in the X-ray spectra of 
active galaxies (AGNs) and quasars, has been notoriously difficult to measure in 
the high-luminosity, radio-loud quasar 3C 273 (z = 0.158). On the few occasions 
that it has been detected its intrinsic width has been thought to be narrow 
(FWHM < 10, 000 km s" 1 ) with an equivalent width (EW) of a few tens of 
eV. This was consistent with the general trend that as one goes from low to 
high luminosity AGNs the Fe Ka line goes from being strong (EW ~ 200 — 300 
eV) and broad (FWHM ~ 10 5 km/s) to being weak and narrow, or absent 
altogether. Here we present the results of new ASCA and RXTE observations, 
together with archival ASCA data, and show for the first time that the Fe Ka 
line in 3C 273 is as broad as that seen in Seyfert 1 galaxies. The line is resolved 
in two of the observations, with a mean Gaussian width of 0.8 ± 0.3 keV, 
corresponding to a FWHM of ~ 0.3 ± 0.1c. The smallest and largest equivalent 
width measured is 43^34 eV and 133lf 3 eV respectively (quasar frame). The 
Compton-reflection continuum is less than 10% of that expected from a centrally 
illuminated semi-infinite, face-on, Compton-thick disk, confirming previous 
studies that Compton reflection is negligible in 3C 273. The largest values of 
the Fe Ka line EW are under-predicted if the line originates in the disk, unless 
a time lag longer than several days between line and continuum and/or an 
over- abundance of Fe is invoked. We cannot unambiguously constrain the disk 
inclination angle. About 60° is preferred for a cold disk, while a face-on disk is 
allowed if the ionization state of Fe is H-like. 

Subject headings: galaxies: active - quasars: emission lines - quasars: individual: 
3C 273 - X-rays: galaxies 
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1. INTRODUCTION 

The Advanced Satellite for Cosmology and Astrophysics (ASCA - Tanaka, Inoue, and 
Holt 1994) found that the Fe Ka fluorescent emission line in Seyfert 1 galaxies is often very 
broad and the data are consistent with an interpretation of an origin from cold matter in 
an accretion disk rotating around a central black hole (see Fabian et al. 2000 and references 
therein). The line profile is sculpted by characteristic gravitational and Doppler energy 
shifts. Currently, study of the Fe Ka emission line is the only way to probe matter within 
a few to tens of gravitational radii of a black hole. These emission features become less 
common in active galactic nuclei (AGNs) with a 2-10 keV intrinsic luminosity exceeding 
~ 10 45 erg s _1 (or synonymously, quasars). On the few occasions that an Fe Ka line is 
detected in high-luminosity AGN it appears to be narrow (FWHM < 10, 000 km/s) with a 
peak energy indicative of highly-ionized matter. This trend is discussed in Nandra et al. 
(1997a,b), Reeves et al. (2000) and George et al. (2000) and references therein. 

The bright, high-luminosity (L(2 - 10 keV) ~ 0.8 - 2 x 10 46 erg s -1 ; 
Ho = 50 km s _1 Mpc -1 and go — throughout) radio-loud quasar 3C 273 (z = 0.158) 
has been extensively studied in all wavebands (e.g. see von Montigny et al. 1997; 
Tiirler et al. 1999 and references therein). In particular, X-ray studies have found a 
featureless nonthermal continuum, with the occasional detection of an Fe Ka emission 
line. Proportional counter experiments such as those on EXOSAT and Ginga could not 
constrain the width of the line and, for a narrow line assumption, found equivalent widths 
of ~ 10 — 50 eV when it was detected (e.g. Turner et al. 1990; Williams et al. 1992). 
Subsequent observations with ASCA again found some weak detections with an equivalent 
width of the order of tens of eV (Yaqoob et al. 1994; Cappi et al. 1998; Reeves et al. 2000). 
Observations with BeppoSAX told a similar story (Grandi et al. 1997; Orr et al. 1998; 
Haardt et al. 1998). 
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In these ASCA and BeppoSAX studies, the Fe Ka line in 3C 273 has been assumed 
or claimed to be narrow. However, on two occasions (1996 July and 1998 June), during 
coordinated ASCA, RXTE and BeppoSAX multi-mission cross-calibration campaigns, we 
found that the Fe Ka line was in fact broad. This is not inconsistent with previously 
published results because when the equivalent width of the line is small it is very difficult 
to constrain the intrinsic width. Results from the 1996 July and subsequent ASCA 
observations have not yet been published. As for the BeppoSAX results (Grandi et al. 1997; 
Haardt et al. 1998), recalibration of the MECS aboard BeppoSAX has recently brought the 
narrow width of the line into question (Grandi, private communication). In this Letter we 
re-examine the old ASCA data, and present results from three new ASCA observations, 
plus two simultaneous RXTE observations, in order to investigate the nature of the elusive 
Fe Ka emission line in this source and we conclude that it is broad, just like that found in 
Seyfert 1 galaxies. 

2. THE DATA AND SPECTRAL FITTING RESULTS 

Between 1993, June and 2000, January, there were a total of six ASCA campaigns 
on 3C 273, resulting in twelve pointed observations. Two observations have simultaneous 
RXTE data (1996 July and 1998 June). Table 1 shows an observation log, exposure times 
and count rates for these observations. We utilize all four instruments aboard ASCA (SIS 
and GIS), fitting the same model to all detectors simultaneously with only the relative 
normalizations independent. All the SIS data were taken in 1-CCD mode. The ASCA data 
were analysed using the methods described in detail in Yaqoob et al. 1998. Note that 
observation B was a calibration observation performed in four parts, placing the source 
sequentially on each of the four SIS chips over a period of ~ 5 days. Observation C was 
a multiple 'snapshot' science campaign carried out over ~ 11 days. Preliminary spectral 
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fitting showed that we could not obtain useful constraints on the Fe Ka line parameters 
from the separate pointings due to insufficient signal-to-noise so we combined the spectra 
from individual pointings for each of the data sets "B" and "C" . For RXTE we used only 
the Proportional Counter Array (PCA) data, which were reduced using methods described 
in Weaver, Krolik, & Pier (1998), except that we used a later version of the spectral 
response matrix (V 7.01) and background model ('L7_240_FAINT'). Only data from PCU 
0,1 and 2, layer 1 were used in this analysis. Note that cross-calibration studies show that in 
the 2-10 keV band PCA fluxes are systematically at least 20% higher than ASCA whereas 
ASCA and BeppoSAX fluxes agree to a few percent Q. We will rely on the ASCA fluxes for 
consistency since we discuss only two RXTE observations. 

Since we are only interested in the Fe Ka line properties for this work, we only utilize 
data in the 1.5-10 keV band for ASCA and 3-10 keV band for the PCA (the latter does not 
have a useful response below 3 keV). First, we fitted each data set with a simple power-law 
plus Galactic absorption (Nh = 1.79 x 10 20 cm -2 ; Dickey & Lockman 1990), but omitting 
data between 4-8 keV (the Fe-K line region). Visual inspection of the data/model residuals 
over the whole energy range then revealed, in most cases, a hump peaking around 5.5 keV 
(observed frame), presumably due to an Fe Ka emission line. The most dramatic of these 
was for observation D, for which the SIS, GIS and PCA data/model ratios are shown in 
Figure 1. It can be seen that the Fe Ka line is very broad, extending from ~ 4 — 7 keV 
similar to the broad Fe Ka lines in Seyfert 1 galaxies. We emphasize that this is not a result 
of instrument calibration problems. The line is broad as measured by five independent 
instruments, out of which three instrument groups (SIS, GIS and PCA) were calibrated 
independently in the relevant energy band. We show below that even the PCA, despite 
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its poor energy resolution resolves the Fe Ka line at 90% confidence. Before discussing 
quantitative measurements of the line parameters we show in Figure 2 similar data/model 
ratios resulting from co-adding all the SISO and SIS1 data in Table 1. This constitutes all 
the SIS data for 3C 273 taken over the entire lifetime of ASCA, amounting to an exposure 
time of ~ 258 ks. The broad Fe Ka line can be seen with very high signal-to-noise. It is not 
meaningful to perform more detailed spectral fitting on the co-added spectra because it is 
known that the spectral index of the continuum is variable over long timescales (Turner et 
al. 1990) so we will not utilize the co-added spectrum further. Note that in Figure 2, at the 
highest energies (above ~ 8 keV), the SIS1 data are systematically higher than SISO. This 
is a known calibration uncertainty which shows up in other sources when the signal-to-noise 
is high []. 

We then fitted each of the eight data sets (six ASCA and two RXTE ) with a simple 
power-law model (including Galactic absorption) and a Gaussian to model the line, utilizing 
all the data in the selected energy bands {ASCA: 1.5-10 keV, PCA: 3-10 keV). There were 
a maximum of three free parameters associated with the line; the line center energy {E c ), 
intrinsic Gaussian width (ope), and intensity (/Fe)- We found that all except observations 
D/ PCA and E/ ASCA had insufficient signal- noise to constrain both the center energy and 
intrinsic width simultaneously so we fixed ope at 0.8 keV (consistent with the best fitting 
values from observations D/ PCA and E/ ASCA ). In cases where the line was particularly 
weak we also had to fix E c (at 6.67 keV, from the best-constrained measurement which was 
from observation E/ASCA ). The spectral fitting results are shown in Table 2, which also 
shows the values of A% 2 compared to the case when no Fe Ka line was included in the 
model. Clearly, the line is detected with a very high level of confidence ( for the two worst 
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cases it is > 98.2% and > 99.975% for observations B and F respectively). Note that in 
observations D/ ASCA and E/PCA the line is actually resolved at the 90% confidence level. 
The best-fitting line energies from ASCA are all ~ 6.6 keV, and the weighted mean, from 
the ASCA observations that could constrain the energy, is 6.62 ± 0.20 keV. The statistical 
uncertainty is still sufficiently large that all ionization states of Fe except H-like are allowed. 

Previous work has shown that a Compton-reflection continuum, common in Seyfert 
1 galaxies, is particularly weak in 3C 273 (e.g. Turner et al. 1990; Grandi et al. 1997). 
Nevertheless, we tested whether Compton-reflection could affect the deduced width of the 
Fe Ka line. Using the ASCA data with the best-determined line width (observation E), 
we added a Compton-reflection continuum to the power-law plus Gaussian model. For 
the most conservative constraints we assumed the primary X-ray continuum illuminates 
a Compton-thick disk viewed face-on. The only additional free parameter is the effective 
covering factor of the disk, R, as a fraction of of the solid angle 2n. The best-fitting value 
of R was 0.0, with a 90% confidence upper limit of 0.20 (four interesting parameters, or 
Ax 2 = 7.779). Thus, the best-fitting line parameters remained unchanged, with the bounds 
on the line width slightly larger, (<jp e = 0.80 ± 0.51 keV). Thus, Compton reflection does 
not affect our conclusion about the broadness of the Fe Ka line. We used the PCA data for 
observation E to obtain even tighter constraints on Compton-reflection. Using the same 
model, and extending the PCA data up to 20 keV, we obtained R = 0.0 with an upper 
limit of 0.094 (90% confidence, A% 2 = 6.251; recall that the line width had to be fixed for 
the PCA/E data). 

As can be seen from Table 2, the smallest and largest measurements of /Fe are for 
observations A and F respectively and do not overlap at the 90% confidence level, yet these 
observations have 2-10 keV fluxes which differ by only ~ 10%. Observations A and F also 
have the smallest and largest EW respectively, which are also mutually exclusive at the 



90% confidence level. However, a clear trend in the line variability is not apparent from 
these data; studies with XMM will clarify this issue. 

Since the broad Fe Ka emission line (Figures 1 and 2) is similar to that seen in many 
Seyfert 1 galaxies (e.g. Nandra et al. 1997a; Reynolds 1997) we fitted the ASCA data 
with a model in which the line photons are emitted in a disk rotating around a central 
Schwarzschild black hole (e.g. Fabian et al. 1989). The parameters are o bs (inclination 
angle of the disk normal relative to the observer), r, (inner radius), r (outer radius), q 
(power- law index characterizing the line emissivity as oc r~ q ), 1^ (intensity), and Ed (line 
energy in the disk rest frame). The fits were insensitive to r« and r so they were fixed at 
6r g and 1000 r g respectively (r g = GM/c 2 ). The data could not constrain q so we fixed 
it at 2.5, a value consistent with the mean found from a sample of Seyfert galaxies (e.g. 
Nandra et al. 1997a). Unfortunately there is too much interplay between Ed and the other 
line parameters so we fixed Ed at 6.4 keV, corresponding to a cold disk. Thus there are at 
most two free parameters in this model of the Fe Ka line (Id and ^obs)- Moreover, we found 
that only observations D and E could constrain 9 b s , so for the other fits we fixed o b s at the 
best-fitting value for observation D. The results are shown in Table 3. Both observations D 
and E give a consistent value for 6 ^ s of 60°, with a lower bound of ~ 45° (the upper value 
is not bounded). This is somewhat surprising because quasars such as 3C 273 are thought 
to be observed face-on (or near face-on). To investigate whether a higher rest-frame energy 
of the Fe Ka line can lower the inclination angle, we repeated the spectral fits with the rest 
energy fixed at the extreme value of 6.97 keV (corresponding to H-like Fe) and obtained 
#obs — 46^27 an d 22^22 degrees for observations D and E respectively. However, the fits 
are slightly worse, with A% 2 between 2-3. Thus we cannot unambiguously constrain the 
inclination angle and require better data, from XMM for example. Note that EW values 
from the disk model (for both values of the rest-frame energy) are virtually the same as 
those obtained from the Gaussian model. It should also be pointed out that in no case 
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does the disk-line model formally provide a better fit than the Gaussian model (but the 
differences in x 2 are n °t significant). 

3. CONCLUSIONS 

From an analysis of the entire database of ASCA observations of 3C 273 and two 
simultaneous RXTE observations, we have shown that the Fe Ka line in this quasar is 
very broad. The line is resolved in two of the observations with a mean Gaussian width of 
0.8 ± 0.3 keV, corresponding to a FWHM of ~ 0.3 ± 0.1c. The 2-10 keV intrinsic luminosity 
ranged from ~ 0.8 — 2 x 10 46 erg s -1 . The smallest and largest equivalent widths we 
measured are 431^ eV and 1331^3 e V respectively (quasar frame) but the data do not allow 
us to identify any clear trend of line variability and its relation to the continuum luminosity. 
The broad Fe Ka line is similar to the relativistic Fe Ka lines commonly found in Seyfert 1 
galaxies, which are thought to be shaped by gravitational and Doppler shifts in the vicinity 
of a black hole. The EW of the lines in Seyfert 1 galaxies is higher however, typically 
~ 200 — 300 eV. The Fe Ka emission line was previously thought to be narrow (FWHM 
< 10,000 km/s) or absent in AGN with a 2-10 keV luminosity exceeding ~ 10 45 erg s -1 . 
Thus, 3C 273 is now the highest luminosity quasar in which a broad Fe Ka line has 
been found. When modeled by a Gaussian, the weighted mean peak energy of the line is 
6.62 ± 0.20 keV, indicative of He-like Fe, although the 90% confidence errors allow all lower 
ionization states. If we model the Fe Ka line in 3C 273 as originating in an accretion disk 
around a black hole, we obtain a disk inclination angle > 44° if the Fe Ka line originates 
in cold matter. If the line-emitting region is highly ionized the disk inclination angle is 
allowed to be smaller, including near face-on for H-like Fe. The disk model also predicts 
that the Fe Ka line should be accompanied by a Compton-reflection continuum, but we 
constrained the magnitude of the latter to be less than 10% of that expected from a face-on, 
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semi-infinite, Compton-thick disk centrally illuminated by the primary X-ray continuum. 
The EW during some of the observations is too large compared to that predicted for the 
observed reflection continuum (e.g. see George & Fabian 1991). Time lags longer than 
several days between line and continuum and/or an over-abundance of Fe could reduce this 
discrepancy. Time-resolved spectroscopy with XMM is needed to make further progress. 

The authors thank the ASCA mission operations team at ISAS, Japan, and all the 
instrument teams for their dedication and hard work in making these ASCA observations 
possible. The authors also thank Paola Grandi for several years of fruitful collaborative 
work on 3C 273 and ASCA/ BeppoSAX cross-calibration, T. Dotani for helpful comments, 
Keith Jahoda for supplying up-to-date response matrices for the PCA, and the RXTE 
GOF for help with data analysis issues. Finally, we thank an anonymous referee for making 
some important suggestions. This research made use of the HEASARC archives at the 
Laboratory for High Energy Astrophysics, NASA/GSFC. 
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Table 1. 3C 273 Observation Log 



Obs a START (UT) END (UT) Exposure (ks) Count Rate (ct/s) 



A/SISO 


08/06/93 19:09 


07/06/93 23:30 


31.3 


B/SISO 








(1) 


15/12/93 12:17 


15/12/93 23:00 


16.4 


(2) 


15/12/93 23:53 


16/12/93 12:10 


19.0 


(3) 


19/12/93 23:45 


20/12/93 10:08 


17.6 


(4) 


20/12/93 10:45 


20/12/93 21:29 


18.2 


TOTAL 


15/12/93 12:17 


20/12/93 21:29 


71.2 


C/SISO 








(1) 


16/12/93 12:11 


16/12/93 19:57 


11.7 


(2) 


20/12/93 22:04 


21/12/93 04:00 


10.2 


(3) 


23/12/93 23:35 


24/12/93 05:05 


10.4 


(4) 


27/12/93 10:37 


27/12/93 17:59 


10.6 


TOTAL 


16/12/93 12:11 


27/12/93 17:59 


42.9 


D/SISO 


17/7/96 20:05 


18/7/96 13:32 


41.9 


D/PCA 6 


17/7/96 23:17 


18/7/96 13:52 


47.3 


E/SISO 


24/06/98 08:30 


26/06/98 10:00 


73.0 


E/PCA 6 


24/06/98 07:45 


26/06/98 08:43 


16.1 


F/SISO 


09/01/00 16:35 


10/01/00 23:32 


45.2 



4.23 

5.44 
5.12 
4.70 
4.89 
4.91 

5.63 
4.76 
4.80 
5.71 
5.23 
2.44 
22.1 
3.67 
35.0 
3.69 



Note. — All count rates are background-subtracted. 
a Although all four ASCA instruments are used in this analysis, we report here only 
the start/stop times, exposure times and 0.5-10 keV count rates for the SIS0 spectra. 
b The RXTE PCA count rates refer to the sum of PCUs 0,1 and 2, layer 1, left and 
right combined, in the 2-10 keV band (see e.g. Weaver et al. 1998 for details). 
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Table 2. Spectral Fits to 3C 273 with a Gaussian Line Model 



Obs r E c (keV) a Fc (keV) 7 Fe EW (eV) F d X 2 / d.o.f. e A X 2 7 



12.4 1643.4/1542 28.4 

16.8 1467.2/1891 13.7 

16.3 1857.9/1757 21.2 

6.9 1522.3/1485 26.3 

8.7 7.8/13 150.6 

10.7 1679.0/1678 27.5 

13.7 17.0/14 25.8 

11.3 1747.0/1641 5.60 



A/ ASCA b 


1 646+ - 110 

l.O4O_ 110 


6 fi9+ - 36 


0.8f 


2 m+ - 82 

z - iu -0.84 


13311 


B/ ASCA a 


i c 7f) +0.007 

i.ofu_ 0007 


6.67f 


0.8f 


i 1O+0.53 
i - i8 -0.52 


5911 


C/ ASCA b 


1 c;oo+0.008 
± - ooz -0.007 


fi c: 4 +0.44 
D - O4 -0.41 


0.8f 


1 61 +0 ' 75 


70+37 
' y -37 


D/ ASCA b 


1 665+ - 012 

l.DDO_ 0012 


6 64+ - 39 
°-°^-0.38 


0.8f 


1 24+ - 51 


1471™ 


D/ PC A c 


1 646+ - 010 


6.35i°;i 


0.741°;^ 


l-36±JiS 


H6l| 


E/ ASCA c 


1 650+ - 008 


6-67±8:S 


o.so±8:t? 


i-oeiL 7 ? 


85_ 47 


E/ PC A b 


1 675+ - 011 
i - D '°-0.012 


6 4S+ - 46 


0.8f 


i n+0.49 
i - ii -0.48 


61 + 27 


F/ ASCA a 


1 R7t;+0-008 
i - D '°-0.009 


6.67f 


0.8f 


60+ - 47 

U.DU_ 048 


43+34 

40 -34 



Note. -- The continuum model used is a simple power law with photon index, T, plus Galactic absorption 
with Nh = 1.79 x 10 20 cm -2 , plus a Gaussian emission line. The emission line parameters are all referred 
to the quasar frame {z = 0.158); these are the intrinsic width, ape, the center energy, E c , the intensity 
(10 photons cm~ 2 s _1 ), Ip c , and the equivalent width (EW). Fixed parameters are denoted by 'f beside 
the fixed value. 

Errors on T are all 90% confidence for one interesting parameter. 

a Errors on the line parameters are 90% confidence for one interesting parameter (A% 2 = 2.706). 
Errors on the line parameters are 90% confidence for two interesting parameters (A% 2 = 4.605). 
c Errors on the line parameters are 90% confidence for three interesting parameters (Ax 2 = 6.251). 
d Flux in the 2-10 keV observed band, in units of 10" 11 erg cm~ 2 s _1 (four-instrument average for ASCA). 
e d.o.f. = degrees of freedom. 

f The decrease in \ 2 when the Gaussian emission line is added to the simple power-law plus Galactic 
absorption model only. 
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Table 3. Spectral Fits to 3C 273 ASCA Data with a Schwarzschild Disk Line Model 



Obs 6 ohs (degrees) I D EW (eV) X 2 / d -°- f - 



A 
B 
C 
D 
E 
F 



59 f 


-i 7 o+0.56 
L -< 3-0.56 


130tg 


1645.4/1543 


59 f 


o 76+ - 46 

u -'°-0.48 


69l£ 


1473.7/1891 


59 f 


1 qq+0.51 
i -°°-0.50 


76±|8 


1859.7/1758 


59±'l 


i-05l8:if 


i4i -68 


1523.5/1485 


ens 


0.951! 


oo+49 
8 ^-32 


1681.2/1679 


59 f 


n p;i+ - 41 

u - di -0.42 


4 o+34 
4z -34 


1747.3/1641 



Note. - Spectral fits to the ASCA 3C 273 data with 
the same continuum model as in Table 1, but replacing the 
Gaussian model of the Fe Ka emission line with a relativistic 
line from a disk around a Schwarzschild black hole (see text). 
The continuum parameters obtained are virtually identical 
to those in Table 1. The disk inclination is 6 ^, s . The line 
intensity (Id) and equivalent width (EW) are both referred 
to the quasar frame (z = 0.158). Errors are 90% confidence 
for one interesting parameter (A% 2 = 2.706) when 9 ^, s is 
fixed (denoted by T) and for two interesting parameters 
(A% 2 = 4.605) otherwise. 
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Figure Captions 

Figure 1 

The Fe Ka line from the 1996 July campaign of 3C 273 (Obs D). The data were fitted 
with a simple power-law plus Galactic absorption, omitting the 4-8 keV data (see text 
for details). Shown are the ratios of data compared to this model for RXTE PCA (top 
inset), ASCA SISO (filled circles) ASCA SIS1 (crosses), ASCA GIS2 and GIS3 (filled circles 
and crosses, respectively, in bottom inset). The line is clearly broad, just like the profiles 
commonly found in Seyfert 1 galaxies. 

Figure 2 

All the data for 3C 273 taken by ASCA over the entire mission (1993-2000) and added 
together, separately for SISO and SIS1. The total exposure time is ~ 258 ks. Shown are 
the data/model ratios when the co-added spectra were fitted with a power law, omitting 
the 4-8 keV band. The broad Fe Ka line is a very prominent feature. Since the statistical 
errors are so small, systematic problems with the SIS+XRT response show up as a hard tail 
in SIS1 relative to SISO. 
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Fig. 1. — See figure caption. 
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Fig. 2. — See figure caption. 



